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Abstract 

The tetrahedron equation is a three-dimensional generalization of the Yang-Baxter equa- 
tion. Its solutions define integrable three-dimensional lattice models of statistical mechanics 
and quantum field theory. Their integrability is not related to the size of the lattice, there- 
fore the same solution of the tetrahedron equation defines different integrable models for 
different finite periodic cubic lattices. Obviously, any such three-dimensional model can be 
viewed as a two-dimensional integrable model on a square lattice, where the additional third 
dimension is treated as an internal degree of freedom. Therefore every solution of the tetra- 
hedron equation provides an infinite sequence of integrable 2d models differing by the size of 
this "hidden third dimension" . In this paper we construct a new solution of the tetrahedron 
equation, which provides in this way the two-dimensional solvable models related to finite- 
dimensional highest weight representations for all quantum affine algebra U q (sl(n)), where 
the rank n coincides with the size of the hidden dimension. These models arc related with an 
anisotropic deformation of the s^(n)-invariant Heisenberg magnets. They were extensively 
studied for a long time, but the hidden 3d structure was hitherto unknown. Our results lead 
to a remarkable exact "rank-size" duality relation for the nested Bethe Ansatz solution for 
these models. Note also, that the above solution of the tetrahedron equation arises in the 
quantization of the "resonant three- wave scattering" model, which is a well-known integrable 
classical system in 2 + 1 dimensions. 
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1 Introduction 



The tetrahedron equation (TE) [1] is the three-dimensional analog of the Yang-Baxter equation. 
It implies the commutativity of layer-to-layer transfer matrices [2] for three-dimensional lattice 
models of statistical mechanics and field theory and, thus, generalizes the most fundamental 
integrability structure of exactly solvable models in two dimensions [3] . 

The first solution of the TE was proposed by Zamolodchikov [1] and subsequently proven 
by Baxter [4]. Later, Baxter [5] exactly calculated the free energy of the corresponding solvable 
three-dimensional model in the limit of an infinite lattice. Next, Bazhanov and Baxter [6] 
generalized this model for an arbitrary number of spin states, N (for the original Zamolodchikov 
model N = 2). The corresponding solutions of the TE was found by Kashaev, Mangazeev, 
Sergeev and Stroganov [7,8]. The other known solutions, previously found by Hietarinta [9] 
and Korepanov [10], were shown to be special cases of [8]. A review of the most recent activity 
related to the TE can be found in the work of von Gehlen, Pakulyak and Sergeev [11]. 

It is worth mentioning that the generalized Zamolodchikov model of ref. [6] also possesses a 
much simpler integrability condition — the "restricted star-triangle relation", introduced in [12]. 
Remarkably, the very same relation serves as the five-term (or the "pentagon") identity for the 
quantum dilogarithm of Faddeev and Kashaev [13] and has been used by Kashaev to formulate 
his famous "hyperbolic volume conjecture" [14,15]. A review of related mathematical identities 
from the point of view of the theory of basic g-hypergeometric series was given by Au-Yang and 
Perk [16] (see also [8,17]). 

As shown in [6] the generalized iV-state Zamolodchikov model has a profound connection to 
the theory of quantum groups [18,19], namely to the cyclic representations of the affine quantum 
group U q (sl(n)) at roots of unity, q N = 1. For the cubic lattice with n layers in one direction 
this model is equivalent to the two-dimensional sZ(n)-generalized chiral Potts model 1 [26, 27] 
(modulo a minor modification of the boundary conditions). In other words, this particular 
three-dimensional model can be viewed as a two-dimensional integrable model, where the third 
dimension of the lattice becomes the rank of the underlying affine quantum group [6] (for related 
discussions see also [28,29]). 

In this paper we argue that the above "three-dimensional structure" is a generic property 
of the affine quantum groups. It is neither restricted to the root of unity case, nor related to 
the properties of cyclic representations. To support this statement we present a new solution 
of the TE, such that the associated three-dimensional models reproduce the two-dimensional 
solvable models related with finite-dimensional highest weight representations for all the affine 
quantum groups U q (sl(n)) with generic values of q. These models were discovered in the earlier 
80's [30-33] and since found numerous applications in integrable systems. They are related with 
an anisotropic deformation of the s/(n)-invariant Heisenberg magnets [34-37]. In the two-layer 
case (n = 2) these models include the most general six-vertex model [38] and all its higher-spin 
descendants. Our approach is briefly outlined below. 

It is well known that the Yang-Baxter equation 

RabRacRbc = RbcRacRab , (1) 

can be regarded is the associativity condition for the L-operator algebra, defined by the "RLL- 
relation" [39] 

L\ A LlfiRab = RabLl^L\, a • (2) 

1 Note, that even the simplest two-layer case (n = 2) corresponds to the famous chiral Potts model [20-23] 
which has attracted much attention over the last twenty years. The most remarkable recent result there is Baxter's 
derivation [24] of the order parameter conjecture of Albertini et al [25]. 
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In the similar way the TE (here it is written in the so-called vertex form) 

RabcRadeRbdfRcef = Reef RbdfRade Rabc , (3) 

can be regarded as the associativity condition 2 for a three-dimensional analog of the above 
"RLL" -relation, 

L\2,a Lis,b ^23,c Rabc = Rabc -^23,c ^13,6 ^12,o • (4) 

Here all operators act in a direct product of six vector spaces V = V\ ®V2 ®V% ® F a ® Fb ® F c , 




Figure 1: Graphical representation of the operators £12,0 an d R a bc- 




Figure 2: Graphical representation of the tetrahedron equation (fljl. The labels "1, 2, 3, a, b, c" 
on the lines indicate the corresponding vector spaces. 

involving the three identical "auxiliary" vector spaces Vi = V, % = 1, 2, 3, and the three identical 
"quantum" spaces F% = J 7 , i = a,b, c. The operator L\2, a acts non-trivially only in V\ ® V2 (g> F a 
and coincides with the identity operator in all the remaining components of V. The other 
operators in (0J) are defined similarly. The graphical representation of Eq.Q is given in Fig^ 
and Fig|5J Suppose next, that T is a representation space of some algebra A. Obviously the 
space JF® 3 is the representation space of the tensor cube of this algebra, A a <8> Ab ® A c , where 
A a = A® id® id, Ab = id®A®id and A c = id® id® A. Then the operator Li2, a , for instance, 

2 This associativity condition is, essentially, the Zamolodchikov's factorization condition for the scattering 
5- matrix of the "straight strings" in the 2 + 1-dimensional space [1]. 
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can be understood as an operator-valued matrix acting in V\ ® V2, whose elements belong to the 
algebra A, 

L\2,a = £12 (v OJ s a ) , (5) 

where v a denotes a set of generators of A a and s a stands for a set of c-valued parameters. 
With this notation Eq.Q can be re- written as 

) Li 3 (v 6 ,s fc ) L 23 (v 

CI *c CI &c ) , (6) 

where the functional map 

KbM = Kbc 4>R~ b \, V^eA®A®A, (7) 
defines an automorphism of the tensor cube of A, 

TZ : A® A® A ^ A® A® A . (8) 

The relation of the form © is sometimes called the tetrahedral Zamolodchikov algebra [10]. 
Since the operators R solve the tetrahedron equation ©, the map solves the functional TE 

n abc {n ade (K bdf {n cef {(/>)))) = n cef (K bdf (n ade {n ahc (</»)))) . (9) 

Further, denote the result of the action of the automorphism (JJJ) on the basis of A <8> A ® A 

as 

= U a bc{Va) , Vj, = U a bc{vb) , v' c = K a b c {v c ) . (10) 

Eq.@ then takes the form of the local Yang-Baxter equation 

^i2(v a ,s a ) L 13 (v b ,s b ) L 23 (v c ,s c ) = L 2 3(Vc,s c ) L 13 (v' b ,s b )Li2(y' a ,s a ) . (11) 

This equation was introduced by Maillet and Nijhoff [40] and further developed by Kashaev, 
Korepanov and Sergeev [41-43]. It is particularly useful for constructing integrable evolution 
systems on a 2 + 1-dimensional lattice and can be viewed as a three-dimensional lattice analog 
of the Lax-Zakharov-Shabat zero-curvature condition in two dimensions (see [44] for further 
details). 

Suppose now, that the algebra A has a "quasi-classical" limit when it degenerates into a 
Poisson algebra V. In this limit all the generators Vj, i = a,b,c, become commutative formal 
variables, and the relation (jllj) becomes a c- number equation which very much resembles the 
usual quantum Yang-Baxter equation (but does not coincide with it). The quasi-classical limit 
of the map (J7J) defines a symplectic transformation of the tensor cube of the Poisson algebra V, 
which, of course, solves the same functional TE Q. 

The approach we employ in this paper can be viewed as a "reverse engineering" of the above 
procedure. We start from the c-number equation Ijlljl with a simple particular Ansatz for the 
matrices Li 2 , L13 and L 23 acting in the product of two-dimensional vector spaces C 2 <8>C 2 ®C 2 . 
Solving this equation we construct a symplectic transformation of the tensor cube of the Poisson 
algebra V, defined by the brackets 

V : {x,y} = 1- xy, {h, y} = y, {h, x} = -x . (12) 

A natural quantization of this algebra lead to the ^-deformed Heisenberg algebra A = TL q (or 
the (/-oscillator algebra [45,46]), 

H q : gyx-g _1 xy = 9-g~ 1 , [h,y]=y, [h,x] = -x. (13) 
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Surprisingly, only minor modifications of the above symplectic transformation and the form of 
the matrices L are required to obtain the quantum map (j7J) for the tensor cube of H q , which 
solves the quantum variant of the equation Finally, we obtain the corresponding solution 

of the TE © by using an explicit form of the quantum map (J7J) to calculate matrix elements of 
the operator R a bc for the Fock representation of the g-oscillator algebra (|T3|) . 

In this paper we briefly present our main results, postponing the details to the future pub- 
lication [44]. The paper is organized as follows. The solution of the local Yang-Baxter equation 
ljll|) with commutative variables (classical case) is given in Sect. 2. The quantum case is consid- 
ered in Sect. 3. A new solution to the tetrahedron equations (jSJ) and (jlj) is presented in Sect. 4. 
The connection of this solution with certain two-dimensional solvable models and the theory 
of quantum group is discussed in Sect. 5. In Sect. 6 we derive the "rank-size" (RS) duality for 
these two-dimensional models. In Conclusion we summarize our results and discuss some open 
questions. 



2 Local Yang- Baxter equation 

Our starting point is the local Yang-Baxter equation where each matrix L depends on two 
pairs of variables, a pair of dynamical variables v = (x, y) and a pair of (complex) parameters 
s = (A,/i). Altogether the equation (fTTj) contains six different sets Vj, v-, i = a,b,c (each L 
depends on its own set v) and three sets of parameters Si, i = a,b,c (they remain the same for 
both sides of the equation) . 

Any matrix acting in the product C 2 <S> C 2 can be conveniently presented as a two by two 
block matrix with two-dimensional blocks where the matrix indices related to the second vector 
space numerate the blocks while the indices of the first space numerate matrix elements inside 
the blocks. With this conversion define the matrix 



L 12 (v,s) 



: 2 ®c 2 , 



(14) 



as follows 



L 12 (v;s) 



/ 1 






where v = (x,y), s = (A,//) and 

k 2 = 1 — xy , k 












Xk 


y 





\fJ,X 


lik 











— A/i 


-» 1 


when 


xy - 



(15) 







(16) 



Now let us consider (jllj) as an equation for v^,v^,v^,, regarding v a ,Vb,v c and s a ,Sb,s c as 
given. It should be stressed that we are seeking for a solution where no constraints on the 
variables Vj, s«, i = a,b, c, in the left hand side of (fTT|) are imposed. This requirement excludes 
the case when Vj = , where the equation (|11|) reduces to the ordinary Yang-Baxter equation for 
the six-vertex free-fermion model 3 (note that the matrix elements of ()15|) satisfy the free-fermion 
condition). 

Detailed considerations show that Eq. ()lH) contains exactly six algebraically independent 
3 In case of the ordinary Yang-Baxter equation the three matrices L12, £13 and L23 are not independent. 
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equations, which can be uniquely solved for six unknowns x[, i = a,b, c, 



x a — T~ ( k c X a T ^aXbUc ) j 2/ a — ^6 "T - ( ^c2/a X a /J, c k a ybX, 



x'b = x aX c + -^—k a k c x h , y' b = VaVc + A Q /i c /c a k c y b , , (17) 

= k, ( k a x c - k c y a Xb ) , y c = ku ( k a y c — Xa^ckcXaVb 

Ha V A a/ U c J Vb\ 

where k' 2 = 1 — a^-y^, in particular, 

kb = k a kjjk c — 2k a k c + k a + /c c - - — X a n c k a k c x a ybX c . (18) 

Elements A:^ and k' c can be found from 

k'aK = k a k h , fc^c = k b k c . (19) 

The signs of k' a , k' b and k' c are fixed by the condition l(TB|) . 

It is not difficult to verify by explicit calculations that the functional operator IjlOjl associated 
with the transformation (|17I18() satisfies the functional TE ©, as it, of course, should. In 
addition, one can verify that the same transformation preserves the Poisson structure of the 
tensor cube of the algebra (|12j) . where the element k = e~~ h / 2 is constrained by the relation (|16|) . 

k 2 = 1 - X y, k 2 = e~ h . (20) 
In other words, the set of Poisson brackets 

{xj,Vj} = 1 -Xjyj , j = a,b,c, (21) 

(where all other brackets are zero) imply the same set of brackets for the "primed" variables 
x j:Vp 3 = a,b,c. 

It should be noted that the general solution of the local Yang-Baxter equation (jll|) with the 
six- vertex- type operators Lij, satisfying the free fermion condition, was found by Korepanov [42]. 
This solution, however, cannot be immediately understood as a transformation of variables 
between the left and right sides of (|llj) , since it contains some spurious degrees of freedom which 
cancel out in each side of this equation separately. We eliminate these "unwanted" degrees of 
freedom and choose a special parametrization (|15|) of the matrices L, such that the resulting 
map (|17ll8j) preserves the Poisson structure of the tensor cube of the algebra (|12)l. The reasons 
why this is possible at all and why does the algebra ()12|) appear here currently remain unclear. 

3 Quantization 

The simplest quantization of the Poisson algebra (|12j) with the constraint (|20j) is the ^-oscillator 
algebra (|13|) with the element k = q h obeying the relation 

k 2 = 1 - yx, k = q h . (22) 
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Li 2 (v;s) 



(23) 



Using the same conventions as in (|15|) . assume the following Ansatz for the operator- valued 
matrices L in 

/ 1 \ 

Ak y 

— q^A/ix //k 

\ -g _1 A// / 

where v = (x, y), s = (A, /i) and the parameter q is same as in (fTT?|) . Obviously, Eq. (|TT|) becomes 
an operator equation. It involves three mutually commuting sets of generators , j = a, b, c, 

satisfying 1)13(1 and 1)22)1 . Just as in the commuting variable case of Sect. 2, one can solve 1)11)1 
for operators x^-, y^, j = a, b, c. Note, that no assumptions on their commutation properties are 
required. The solution is: 



A a /A 



-k a x b y ( 



y« = k 6 ^ I Ky a — — Kyb*< 



where 



A a /A 



■k a k c x b , 



k, ( k„x^ 

Ha 



AoM. 



y D = y a y c + A a /4ck k c y& , 

/ i 1—1 Ha ( , ^aHc 

yc = K —\ k ayc-—— 

Hb\ Q 



K* a yb 



(24) 



kl 2 



7 2 k 2 k 2 k 2 - (1 + g 2 )k 2 k 2 + k 2 + k 2 - kak f QXfcyC " A a /A,k a k c x a y b x c (25) 

AaHc 



and 



k' k' — k W k' — k k 



(26) 



Again one can verify directly, that this operator transformation is an automorphism of the tensor 
cube of (/-oscillator algebra 1)13122)1 . and that the functional operator ()1(J)) corresponding to this 
transformation solves the quantum functional TE Q as well. 



4 A solution of the tetrahedron equation 

Consider the Fock representation, J 7 , of the g-oscillator algebra ((T3*() spanned on the basis \a), 
a = 0, 1, 2, . . . , oo, 

T : x|0)=0, y\a) = (1 - q 2a+2 )\a + 1) , h\a)=a\a), (a\a) = 5 a ^ (27) 

and assume that q is not a root of unity. Then the representation 1)27)1 is irreducible and therefore 
the equation determines the matrix elements of the operator R a bc to within an unessential 
scale factor. Note, that the constraint (j2Q|) for the representation ((77)1 is fulfilled automatically. 
Explicitly, Eq.0 reads 

Xj- Rabc = Rabc Xj , y^ i?a6c = -Rabc Yj , j = a, 6, C , (28) 

where x^-,y'- are given by (|24(1. Solving these equations one obtains, 

(a, / 9, 7 |i2 a6c |a', / (3',7 / ) = (A C /A 6 ) Q ' (A^c/ W/i 6 ) 7 ' (a, /?, 7 | r |a', /?', 7') (29) 
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where matrix indices (a, a'), (/?,/?') and (7,7')! refer to the representation spaces T a , Tb and 
F c , respectively, and take the values 0, 1, 2, . . . , 00. Note that the parameters Sj = (Xj,fij), 
j = a,b,c, associated with these spaces in 1)1 enter (|29[) only through simple power factors. 
The constant matrix (. . . | r | . . .) is given by 

(a'-/3)( 7 '-/3)-/3 

(a,/?,7|r|a',/?',7'} = R^+^+^+^+y - . 2 2 , Pp<{q 2a , q W , q 21 ) , (30) 

\Q >9 )p 

n-l 

where (x;q 2 ) n = TJ(1 — q 2 ^x). The polynomials P n (x,y,z) are defined recursively, 

P n (x, y, z) = (l- x)(l - z)P n ^(q- 2 x, y, q~ 2 z) - q 2 - 2n xz(l - y)P n - 1 (x, q- 2 y, z) , (31) 

with the initial condition Po(x,y,z) = 1. The first non-trivial polynomial reads 

Pi(x,y,z) = 1 -x -z + xyz , (32) 

whereas all higher polynomials soon become very complicated and their explicit form is not very 
illuminating 4 . 

The recurrence (|3*Tj) ensures that all the defining relations (|28j) for the operators R a bc (which 
are equivalent to the TE (jlj)) are satisfied. Thus the operators L and R a bc defined by (|5|) . (f23|) 
and Q29|) . respectively, solve the TE @. Further, the general arguments of Sect.l and the 
irreducibility of the representation (|27|) imply that the operators R a bc, defined above, also solve 
the TE ©. 

Note, that all the parameters Xj and fij associated with six different vector spaces labelled by 
a, b, c, d, e and / totally cancel out from Eq.JSJ). Therefore, this equation is essentially a constant 
TE equation for the parameter- independent matrix r, ()3U|) . 

r a fec Tade *bdf ?cef = ?cef *bdf ^ade *abc , (33) 

which coincides with Eq.(|SJ) when all the A's and /i's therein are equal to one. The same 
remark applies to the other TE Q (and also to Ea. (|34|) below) as well. Nevertheless, the above 
parameters play an important role for the associated two-dimensional i?-matrices considered in 
the next section and that is why they are retained here. Note also, that due to the presence of two 
delta functions in (|30|) each side of (|33|) contains only a finite summation over the intermediate 
matrix indices. 

There exists another relevant TE, which we present without a proof. It involves by the same 
operators L, 

M o0 , a (/V//) M lU (A'/A) L 01>b (s)L- oU (s') = L m ^s')L Qhb ( s )M 1IjB (A'/A) M^ a (/*//*') (34) 

and new operators Mj Jja which are very much similar to L's. They act non-trivially in the 
product of three spaces Vi Vj ® T a , where V { = Vj = C 2 and T a is the Fock space ((23) • With 
the same conventions as in (|15p these operators are defined as 

/ A h \ 

kA h yA h 

g- x xA h kA h 

\ q- x \^ J 

where x, y, h are the generators of g-oscillator algebra (|13|) (recall that the element h has the 
integer- valued spectrum (|27|). so the power A h is well defined). The proof of (|3*4"|) is given in [44]. 

Nevertheless, V.V.Mangazeev [47] found an elegant closed form solution for these polynomials in terms of the 
hypergeometric function, which will be reported elsewhere [48]. 



Mi2,«(A) 



(-q) h , (35) 
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5 Quantum R-matrices for U q (sl(n)) 

Consider a model of statistical mechanics on a cubic lattice with the toroidal boundary conditions 
in all three lattice directions. Let each edge carry a discrete spin variable taking the values 
0.1.2.... , oo and each vertex is assigned with Boltzmann weights, given by the matrix R a b c , 
(|29|) , so that its indices are identified with the spin variables on six edges surrounding the vertex. 
Namely, the three pairs of indices (a, a'), /?') and (7, 7') are associated with the edges oriented 
along three different lattice directions "a" , "6" and "c" respectively. 




Figure 3: A chain of the operators R traced in the space T t 



It is well known that any edge-spin model on the cubic model can be viewed as a two- 
dimensional model on the square lattice with an enlarged space of states for the edge spins 
(see [2,6] for additional explanations). Consider the following quantity, 

Rbc = Trjr (R a biciRab2C2 ■ ■ ■ Rab n c^\ , (36) 

which involves the product of the vertex weights along a whole line of vertices in the direction 
"a" of the lattice, with n being the corresponding lattice size (see Fig. [5J). It is an operator 
acting in the product of two "composite" spaces Fb ® F c , 

F b = :F 6l ®.7 r 6 2 ® "-®:F& n , F c = T cx <g> F C2 <8> • • • ® T Cn , (37) 

where each space T coincides with the Fock space (|27[). It is easy to show that the operator 
(|3*6*|) satisfies the Yang-Baxter equation, 

RbcRbdR-cd = R-cdR-bdR-bc (38) 

as a simple consequence of the TE © and the fact that the matrix ()3U|) is a non-degenerate 
matrix in .F® 3 '. Evidently, the operator 1)36(1 is the i?-matrix of the associated two-dimensional 
model. Similarly, one can construct the L-operator 

Lvb = Try 1 £01,61-^02,62 ' " " Lon,b n ) • (39) 
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where L^j, is given by ((15)1 and trace is taken over the two-dimensional space Vq = C 2 . Eq.@ 
implies the Yang-Baxter equation 

LvbLvcR-bc = R-bcLvcLvb • (40) 

The local 3d structure of this equation is shown in Fig. 0J Evidently, the i?-matrix (|36l) 
"intertwines" the L-operators (|39|) in the quantum spaces. There exists another i?-matrix, 
Rvv') which intertwines the same operators in the auxiliary spaces 

LvbLv'bRvv = Rvv'Lv'bLvb • (41) 
This Yang-Baxter equation follows from (|34() . 




Figure 4: The local 3d structure of the product LvbLvcR-bc m the left hand side of the Yang- 
Baxter equation (|4*0]> . 

In general, each operator Loj,b 3 m P^j) depends on its own set of parameters Sj = (Xj,/j,j). 
Consider the case when all Xj are the same, Xj = X, j = 1, 2, . . . , n, whereas the parameters 
{/j,} = {fj,i, /J,2, • • • > /%} are kept arbitrary; we will write the operator (|39|) as Lvb(A|{/u}) to 
indicate these arguments explicitly. This operator acts in product V (g) Ft,, where 

V = (C 2 )® n (42) 

and Fb is defined in (|3*7j) . It can be regarded as a 2 n by 2™ matrix with operator- valued entries 
acting in the quantum space Fb- It turns out, that this matrix has a block-diagonal structure, 
with (n + 1) blocks of the dimensions, which coincide with dimensions of the fundamental 
representations ■K UJk of the algebra U q (sl(n)), 

dim(7r^) = k] ( n i k y , k = 0,l,...,n. (43) 

Here ut denotes the highest weight of the A:-th fundamental representation. Note, that ir^o an d 
iTu n are the trivial one-dimensional representations. Consider, for instance the n-dimensional 
subspace of (|42|) with the basis 

m = (j)®---(o)®(l) <8 (o)®"' (o)' i = ^---^. (44) 

i th place 
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One can easily see that the operator (|39|) acts invariantly in this subspace. The corresponding 
n by n matrix reads 

(^|Lvb(A|M)|^'} = Ai(A|{/4) = -q^mbjW, (45) 
where Cij(X) does not depend on the parameters 

n 

£(A) = Y J E aa ® (X n q J ~ ha - q h -) + Y^ (x^ a E a p <g> Sp a + \ n+a ^Ep a ® q J 8 a p) . (46) 

a=l a</3 

Here (E a p)ij = S^S^p is the standard matrix unit and 

n 13 13-1 

7=1 7=a 7=q+1 

where (x a , y a , h a ), a = 1,2, ... , n are the generators of the n independent copies of the q- 
oscillator algebra (|13|) . A very simple inspection shows that the matrix C(X), given by (|46[) . is 
nothing but the fundamental L-operator [19,49] for the quantum affine Lie algebra U q {sl{n)) 
where the parameter A is just the usual spectral parameter. Remind, that the spectral parameter 
dependent L-operators arise as specializations of the universal i?-matrix [18, 50, 51] when the 
affine algebra U q (sl(n)) is realized by means of the evaluation homomorphism [19] into the 
finite-dimensional algebra U q (sl(n)). A specific feature of the L-operator ()46j) is that the latter 
algebra is realized in this case with the help of another homomorphism into the tensor power of 
the g-oscillator algebra (|T3)) . 

The Cartan-Weyl generators E a , F a , H a , of the algebra U q (sl(n)), where a = 1, 2, ... n — 1, 
satisfy the commutation relations 

[H a ,Ep] = A a> pEp , [H a ,Fp] = —A aj pFp , [E a ,Fp] = 5 aj p — - — — (48) 

and the cubic Serre relations 

E a E a±l - (q + q' 1 ) E a E a±1 E a + E a±1 E 2 a = , 

F a F a±l ~ (4 + Q' 1 ) F a F a±l F a + F a±l F a = > 



(49) 



where A a p = 25 at p—5 a+ i t p—5 a: p + i is the Cartan matrix. The L-operator corresponds to the 
following well-known realization of these defining relations in terms of the (/-oscillators [45,46], 

1 h a + l ha 1 _h 2L h Q + l 

Ea = ( g _ g -l) 9 2 S <*,*+ iq 2 = ( q -q-i) Q 2 y ^a+iQ 2 , 



1 h q + l _ha I h a + l hg (§()) 

Fa = (q-q-±) Q 2 £a+1 > aq 2 = ( q -n-l\ q 2 y<x+l x a Q " 



H a = h a - h a+ i , a = 1, . . . , n - 1 . 

This correspondence immediately follows from a comparison of Q46JI with the general form of 
the fundamental L-operator for the algebra U q (sl(n)) given in [19,49]. 

Note, that the parameters {fi} enter Ea. (|45|) trivially; their only effect reduces to the pre- 
multiplication of £(A) by the diagonal matrix diag(/zi, /X2, . . . , fj, n ). This corresponds to the 
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introduction of arbitrary "horizontal fields" [33], which, obviously, does not affect the integra- 
bility of the 2d model related with the L-operator (|46|). 
The general decomposition of (|39|) has the form 

n 

Lvb(A|M) = C sl ^(oo k ,X\M), (51) 

k=0 

where C sl( - n \oJk, A |{/i}) is the L-operator corresponding to the fc-th fundamental representation 
ir Wk of U q (sl(n)) (with the horizontal fields {//}); in particular, £ sl( - n \uJi, A |{/i}) exactly coincides 
with C(X\{fik}) defined by (|45J46|) , For further references define row-to-row transfer matrices 
corresponding to the above L-operators for an inhomogeneous periodic chain of the length m, 

TW(u k , {A} |{ M }) = Tr^ (c'Wfa, A m |{ M }) C sl ^(u; k , A m _! |{/x}) • • • £ S ^W Aj |{m})) , 

(52) 

where {A} = {Ai, A2, . . . , A m }, denotes the set of spectral parameters. 

For the Fock representation ()27|) the element JT, defined in (|47() . has the integer- valued 
spectrum J = 0, 1, 2, . . . , 00, where J is just the total occupation number of n (/-oscillators. 
Thus, the quantum space T® n can be split into a direct sum 

F*n = ® Vj , dimVj= ( J + n -^ (53) 
j=o J\(n — iy. 

where the spaces Vj are isomorphic to the representation spaces ttj UJi of the algebra U q {sl{n)) 
(the J-th symmetric powers of the vector representation 7r wi ). As the L-operator (|3H|) commute 
with the element J , it acts invariantly in each component of Indeed, the L-operators ap- 
pearing on the right hand side of (|5l"j) split into a direct sum of U q {sl{n)) i?-matrices, Hjui,u h , 
corresponding to all representation ir j Ul , J = 0, 1, 2, . . . , 00 in quantum space. Likewise, the op- 
erator (|36|) splits into an infinite direct sum of (appropriately normalized) i?-matrices, TZj j 1 j'uj 1 , 
corresponding to the representations ir j Ul and irj/ Ul 

00 00 

(54) 

J=0 J'=0 

where, for brevity, the spectral parameter and field arguments are suppressed. This decompo- 
sition can be established directly from the definition (|36|k the calculation of each term in the 
right hand side requires the use of the polynomials P n (x,y, z), recursively defined by (|31[). up 
to the order n = J (see [44] for further details). 



6 The rank-size duality 

Consider now the layer-to-layer transfer matrix, T mn of m columns and n rows obtained as a 
trace of a product of nm operators (J23J), 



T mn ({A}|M) = Trace 



(55) 



where the ordered products are defined as 



fj = f /1/2 ■■■ fn , Y[ 9i 9mg m -i "'91 . 



def 



(56) 
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Figure 5: The spacial structure of the layer-to-layer transfer matrix (|55|). 



The spacial structure of T mn is illustrated in Fig[5J Each node located at the intersection 
of the i-th column and j-th row corresponds to the operator Lij(vij, Sij) in (|55|). as shown in 
FigEJ This operator acts in its own quantum space Tij (these spaces are associated with the 
vertical lines in Fig|SJ) and in two two-dimensional vector spaces C 2 which are associated with 
the column and row corresponding to this node. The trace in (|55j) is taken over all these two- 
dimensional spaces. Further, each column is assigned with the parameter Aj, i = 1, 2, . . . , m and 
each row is assigned with the parameter fj,j, j = 1,2, ... , n. The arguments S{j used in (|55|) are 
given by s„ = (A*,/^). 



Liji^ij , Sij) — 3 



Figure 6: Graphical representation of Ly (vjj , s ij ) in 

The layer-to-layer transfer matrices T mn ({A} act in the direct product J : ® mn of the 

Fock spaces (|27|) . They form a two-parameter commutative family of operators. Indeed, the 
Ea.(|34|) implies that 

[T m „({A} T m „({uA} \{vn})] = 0, (57) 

where 

{uX} = {u\ 1 ,u\ 2 , ■ ■ ■ ,u\ m }, {vfi} = {vm,vn2, ■ ■ ■ ,v/jL n }, (58) 
with u and v being arbitrary complex parameters. By construction, the layer-to-layer transfer 
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matrix 

m n 

T m „({uA} \{ Vf x}) = G ^ ( 59 ) 

k=o e=o 

is a polynomial in u and v. Its coefficients, Gke, provide a set of mutually commuting "integrals 
of motion" for this 3d model. The latter non-trivially depend only on the ratios of A's and the 
ratios of /x's, since their overall normalizations can be absorbed into the parameters u and v. 

The order of the products in the definition l)55|) can be interchanged by using a simple 
symmetry property, 

Li2(v,(A,/i)) = [L 21 (v,(/i,A))]* , (60) 

where the superscript u t" denotes the transposition of the four by four matrix (|23j) . In this way 
one obtains the following "reflection" symmetry relation 

T mn ({A} |M) = T nm ({/4 |{A}) . (61) 

Substituting now the definition (|39j) for the inner products in two variants of Eq, (|55j) (for each 
side of the last equation), and using the decomposition (|51|) and the definitions ()52j) . one arrives 
at a remarkable duality relation for the transfer matrices of the associated two-dimensional 
models 

n m 

£ v km , {u\} IM) = £ u fa Tf m \u t , { W } |{A}) . (62) 

k=o e=o 
The left and right hand sides of this relation are just different expansions of the same layer- 
to-layer transfer matrix (|59|) . It provides two different representations for the same integrals of 
motion G^i m terms of two-dimensional lattice models related with two different quantum affine 
algebras U q (sl{n)) and U q {sl{m)). The above relation simultaneously interchanges the rank of 
the quantum algebra and the set of horizontal fields in one model with the length of the chain 
and the set of (inhomogeneous) spectral parameters in the other. We will call this relation the 
"rank-size" (RS) duality relation. Note, that q — > 1 limit of this relation covers the 2d solvable 
models with the Yangian symmetry [18] (the s/(ra)-invariant Heisenberg magnets and associated 
models [34-37]). 

The zero-field case (i.e., when all A's and all /x's are equal to one) corresponds to the ho- 
mogeneous models in both sides of the relation ()62|) . with the spectral parameters u and v 
respectively. For example, when n = 2 it relates the (zero-field homogeneous) 6-vertex models 
on the chain of the length m with an U q (sl(m))-type model [30-33] on the chain of just two 
lattice sites. 

It is well known that all two-dimensional models discussed here can be solved by the nested 
Bethe Ansatz [35-37,52], which allows one to obtain the exact expressions for the eigenvalues 
and eigenvectors of the transfer matrices in terms of solutions of certain algebraic equations, 
called the Bethe Ansatz equations. The above RS-duality relation implies a complete equiva- 
lence between two different nested Bethe Ansdtze associated with the two sides of (|62j) . This 
remarkable (and rather unexpected) equivalence certainly deserves further studies, which are 
postponed to the future publication [48]. 

7 Conclusion 

In this paper we have constructed a new solution for the set of the interrelated TE's ©, Q and 
(|34j) . These equations involve the two different spaces of states for the discrete edge "spins": the 
two-dimensional vector space C 2 and the infinite-dimensional Fock space, J 7 , defined by (|27|) . 
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In particular, the 3d L-operators Ljj ja , entering Eq.@, act in the direct product of the two 
spaces C 2 and one Fock space J- a . The most important ingredient of our approach is the local 
Yang-Baxter equation 1)11(1 with the specific Ansatz 1)23(1 for these 3d L-operators. This equation 
defines the automorphism ()24J) of the tensor cube of the g-deformed Heisenberg algebra (|13|) . 
which almost straightforwardly provides the solution the TE's ((3J) and Q. 

Currently, it is unclear whether the same scheme should work with different 3d L-operators, 
other than 1)23)1 but it is certainly worth exploring this point. Our choice of the form of this 
operator was motivated by a consideration of the "resonant three- wave scattering" model [53,54], 
which is a well known example of integrable classical system in 2 + 1-dimensions. This model is 
described by the non-linear differential equations, 

d a Aj3j(x) = A/3 a (x)A ai (x), (a,/3,7) = perm(l,2,3) , (63) 

for the fields A a n{x), a ^ (3, in the three-dimensional space with the coordinates x = (x\, X2, £3) 
and d a = d/dx a , a = 1, 2, 3. They can be thought as the consistency conditions for the following 
auxiliary linear problem, 

daVpix) = A af3 (x) * a (z), a ^ 13 . (64) 

The 3d L-operator (|23|) naturally arises in the quantization [44] of a discrete lattice analog [42] 
of this linear problem. 

As explained in Sect. 5, every solution of the TE provides an infinite series of quantum R- 
matrices for 2d solvable lattice models (i.e., the solutions of the Yang-Baxter equation). Our 
solution of the TE reproduces in this way the .R-matrices related with the finite dimensional 
highest weight evaluation representations for all the quantum affine algebras U q (sl(n)) with 
n = 2,3,... ,00. This 3d interpretation provides a completely new insight into the properties 
of 2d solvable lattice models, in particular, it leads to a remarkable rank-size (RS) duality, 
discussed in Sect. 6. 

Plausibly, a similar 3d interpretation also exists for the trigonometric i?-matrices related 
with all other infinite series of quantum affine algebras [55,56] and superalgebras [57]. If so, the 
corresponding 3d models will, most likely, involve non-trivial boundary conditions [58-60]. We 
hope to address this interesting problem in the future. 
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